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ABSTRACT 

About 500d after explosion the light curve of the Type la SN 1998bu suddenly flattened and at the 
same time the spectrum changed from the typical nebular emission to a blue continuum with broad 
absorption and emission features reminiscent of the SN spectrum at early phases. We show that in 
analogy to SN 1991T (Schmidt et al. 1994), this can be explained by the emergence of a light echo 
from a foreground dust cloud. Based on a simple model we argue that the amount of dust required 
can consistently explain the extinction which has been estimated by completely independent methods. 
Because of the similar echo luminosity but much higher optical depth of the dust in SN 1998bu compared 
with SN 1991T, we expect that the echo ring size of SN 1998bu grows faster than in SN 1991T. HST 
observations have indeed confirmed this prediction. 

Subject headings: supernovae: general — supernovae:individual(1991T,1998bu) — dust — extinction 



1. INTRODUCTION 

The Type la SN 1998bu was discovered on May 9.9 
(N ; UT by Villi (1998) in a spiral arm of NGC 3368 (M96), 
^J" . a nearby Sab galaxy. It had very good photometric 
C*~) 1 and spectroscopic coverage by several groups (Jha et al. 
1999; Hernandez et al. 2000; Suntzeff et al. 1999) show- 
ing that the luminosity decline was slower than the av- 
erage (Ami5(B), the magnitude decline in the early 15 
days after maximum, was 1.01 ± 0.05 mag). This is al- 
■ most as slow as the slowest SN la on record, SN 1991T, 
O-i] although, unlike SN 1991T, SN 1998bu was spectroscopi- 
q . cally normal before and around maximum. The distance 
!— i ' to the host galaxy was measured using Cepheids variables 
^ i (Tanvir, Ferguson & Shanks 1999) as \x = 30.25 ± 0.19. 
d • However, a calibration of the absolute luminosity of the 
^ ' SN depends also on the estimated reddening which in the 
• i— i , case of SN 1998bu is significant. Hernandez et al. (2000) 
^ ■ estimated E(B -V) = 0.86 ± 0.15 based on the compari- 
J-j \ son of the spectral energy distribution of SN 1998bu with 
. *P_ ! that of the template SN la 1981B. 

SN 1998bu was included in our program of monitoring of 
the nebular phase of SN la. During this monitoring we re- 
alized that SN 1998bu was not declining with the expected 
rate; this, combined with the change in the spectral ap- 
pearance gave the first evidence of the emergence of a light 
echo in SN 1998bu (Cappellaro et al. 2000). These obser- 
vations are presented here for the first time along with a 
simple modeling. 

2. EVIDENCE FOR A LIGHT ECHO IN SN 1998BU 



When the SN la ejecta become optically thin, about 100 
days after the explosion, the SN luminosity is determined 
essentially by the deposition of the kinetic energy of the 
positrons released by the decay of 56 Co. There is evidence 
that even positron deposition may not be complete (Cap- 
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pellaro et al. 1997; Milne, The & Leising 1999), and so 
at advanced phases the light curve declines either at the 
56 Co rate (0.98 mag/lOOd) or faster. The only previously 
known case of a SN la with a decline slower than the 56 Co 
rate was SN 1991T, whose slow late decline and peculiar 
spectra were interpreted as due to a light echo (Schmidt 
et al. 1994). 

A first indication that SN 1998bu was beginning to de- 
viate from the normal exponential decline came on De- 
cember 4 1999, roughly 500 days after maximum. At this 
epoch the observed magnitude (V = 20.75 ± 0.08) was al- 
ready 2 mag brighter than the extrapolation of the radioac- 
tive decay tail. Additional photometric observations over 
the following few months showed that the SN luminosity 
remained almost constant (the decline was only ~ 0.5 mag 
in 100 days). Also, the color at these epochs was unusually 
blue (B — V ~ —0.1). This behavior was reminiscent of 
that exhibited by SN 1991T (Schmidt et al. 1994; Sparks 
et al. 1999). 

Fig. 1 shows the absolute V light curves of SN 1998bu 
and SN 1991T (Cappellaro et al. 1997). Both light curves 
have been calibrated using the Cepheid-based distances 
(Saha et al. 2000; Tanvir, Ferguson & Shanks 1999) and 
corrected for total extinction (cf. Tab. 1). The strik- 
ing similarity of the two light curves suggests that the 
same mechanism causing the flattening of the light curve 
of SN 1991T was also acting in the case of SN 1998bu. 

Although this is not clearly visible in Fig. 1, the abso- 
lute magnitude of SN 1998bu at maximum was about 0.5 
mag fainter than that of SN 1991T. Based on the Ami5 (B) 
vs My relation (Phillips et al. 1999) the difference should 
be only 0.05 mag. The gap increases slightly with time, 
and it reaches about 0.8-0.9 mag one year after maximum. 
That SNe la with a very similar Am^S) actually show 
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significant dispersion in some of their other properties is 
a well-known fact (Mazzali et al. 1998), and is one of the 
main caveats for the use of SN la as distance indicators. 
However, this is not essential for the discussion presented 
here. 
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FIG. 1. The light curve of SN 1998bu is compared with that of 
SN 1991T (circles are our own measurements, whereas triangles are 
data from literature). For reference, the continuous line represents 
the V absolute light curve of SN 1996X, a normally declining SN la 
(the star marks the epoch of the last measurement of SN 1996X and 
after that the line is simply an extrapolation). The dotted line is the 
fit of the observations with the simple model described in Sect. 2. 

The slow photometric evolution of SN 1998bu makes it 
a viable target for spectroscopy even 2 yrs after explosion. 
Spectra of SN 1998bu where obtained using the ESO 3.6m 
telescope (and EFOSC2) at La Silla at two epochs, 2000 
March 3 and March 30. We used grism #11 which, com- 
bined with the 1.2 arcsec slit, allowing one to cover the 
range 3300-7500 A with a resolution of about 15A. Even- 
tually, since no spectral evolution was apparent between 
the two epochs, the two spectra were merged to improve 
the S/N. The total integration time was 9000 sec. The re- 
sulting spectrum, labeled with the average phase of 670d, 
is shown in Fig. 2. For comparison we also show the spec- 
trum of SN 1991T obtained with the same telescope at a 
similar phase. We also show spectra of the two SNe at a 
much earlier epoch, when the luminosity decline was still 
tracking the radioactive exponential tail. 

The similarity of the spectral evolution of the two SNe is 
remarkable. About 300 days after maximum the spectra 
are dominated by the strong emission lines of Fell] and 
Felll] which are typical of SNe la in the nebular phase. 
These emission lines originate in the iron nebula and are 
powered by the radioactive decay of 56 Co to 56 Fe (Kuchner 
et al. 1994). 

In correspondence to the flattening of the light curve 
the spectra change completely, showing a blue continuum 
with superposed broad absorption and emission features. 
These spectra resemble the photospheric epoch spectra, 
but they do not match any specific early-time epoch. For 
SN 1991T it was convincingly suggested (Schmidt et al. 



1994) that the peculiar late time light curve and spec- 
trum can be explained assuming that the emission at these 
epochs is dominated by the echo from circumstellar dust 
of the light emitted by the SN near maximum. Our obser- 
vations suggest that the same mechanism is at work also 
in SN 1998bu. 

The key clement for such a mechanism to work is the 
presence of a dust layer in the vicinity of the SN, scatter- 
ing towards the observer a fraction of the light emitted by 
the SN. Because of the different travel times the observed 
spectrum is a combination of the early time scattered spec- 
tra and of the direct, late-time spectrum. In principle, the 
details of the echo (intensity, duration, spectrum) depend 
on the spatial distribution of the dust and on the physics 
of the scattering process. A careful analysis and modeling 
of time distributed observations may allow one to derive 
useful constraints. Such a detailed study is the subject of 
future work (Patat et al. in preparation), while here we 
illustrate the basic principles of the phenomenon. 

In general, the light curve of the echo is related to the 
light emitted by the SN by the following relation (Cheva- 
lier 1986): 



F echo (t) = f F SN (t - t') f{t') dt' 
Jo 



(1) 



where FsN(t — t') is the flux of the SN at time t — t' 
and f(t), the fraction of this light which is scattered to- 
wards the observer, depends on the geometry of the sys- 
tem and on the properties of the dust (units of f(t) are 
s ; cf. Equation 2). As a first approximation we can 
assume that the SN light curve is a short pulse with a 
duration AtsN, during which the emitted flux F$n is 
constant. Under this assumption Equation 1 becomes 
Fecho(t) = Fsn f(t) AtsN, where At$N can be obtained 
by a numerical integration of the observed light curve as 

Fsn^sn = Jo + °° FsN(t)dt. In the case of the V light 
curve of SN 1991T (Schmidt et al. 1994; Sparks et al. 1999) 
this gives AtsN~0.11 yr. 

As one of the possible geometries which are consistent 
with the observations, we consider a thin dust sheet lying 
in front of the SN and extended perpendicularly to the 
line of sight. For the sake of simplicity we make the hy- 
pothesis that the sheet thickness AD is much smaller than 
the distance between the SN and the dust elements and 
we consider only single scattering. Given this idealized 
configuration and assuming the scattering phase function 
by Henyey & Greenstein (1941), an analytic expression for 
the function f(t) can be derived relatively easily (Cheva- 
lier 1986; Xu, Crotts & Kunkel 1994) as: ' 



1-.9 2 



8^ D + Ct (l +ff 2_ 2^^)3/2 



(2) 



where D is the distance of the dust layer from the SN, Td 
is the optical depth of the dust layer, Lo d is the dust albedo, 
and g is a parameter which represents the degree of for- 
ward scattering. The value of g ranges from for isotropic 
scattering to 1 for purely forward scattering. Empirical 
estimates and numerical calculations (White 1979) give 
g «0.6. For the dust albedo we have adopted u>d ~0.6 
(Mathis, Rumpl & Nordsiek 1977). 
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Comparison of the spectra of SN1998bu (thick lines) and SN1991T 
(thin lines) at two epochs, roughly 1 and 2 years past explosion. 
Spectra have been reported in the host galaxy rest frames, but have 
not been corrected for reddening 

Equation 2 shows that if the dust is too close to the 
object, even if it survives the exposure to the enormous 
radiation flux of the SN, it would not produce a constant, 
long duration echo. This is because of both the angular 
dependence of the scattering phase function and the SN 
radiation dilution. In the case of SN 1991T the echo lu- 
minosity declined by a factor of 4 over ~ 7 yrs (Sparks 
et al. 1999) which implies D > 45 lyr (equation 2). On 
the other hand, if a dust layer is to produce a significant 
light echo, it must be located close enough to the SN that 
the emission from the SN is not too diluted. In fact, other 
conditions being similar, for a distant cloud (D 3> ct), the 
echo flux is F ec ho oc t^/D. In particular, since we found 
that the cloud must be located farther than 45 lyr from the 
SN, we derive a lower limit for the dust optical depth of 
Td > 0.28. Since in our model the dust layer is located in 
front of the SN we expect that it causes some extinction 
on the line of sight. Assuming that the column density 
is constant along different lines of sight and based on the 
relation Ay — 1.0867V (Cox 2000) we can estimate for 
SN 1991T a lower limit for A v > 0.3 mag. 

Actually, both SN 1991T and SN 1998bu show evidence 
of significant extinction from interstellar dust associated 
with the host galaxy, A^ ost = 0.46 for SN 1991T and 
Ay ost = 0.86 for SN 1998bu (Tab. 1). These values were 
derived mainly by photometric methods based on the anal- 
ysis of the light curve and the SN color, and therefore they 
are affected by large uncertainties. In a simplified scenario, 
one can assume that there is only one cloud which is caus- 
ing both the observed extinction and the light echo. In 



this case, since we fix the column density of the cloud, we 
can derive the distance D of the cloud itself from the SN 
through a fit to the echo luminosity. The light curve fit 
for SN 1991T, shown in Fig. 1, gives D ~ 120 lyr. In 
the case of SN 1998bu the light curve is similar to that 
of SN 1991T, hence At SN ^ 0.1 yr, but the optical depth 
of the dust is about a factor 2 larger: therefore we obtain 
D ~ 230 lyr. Note that these estimates are independent 
of the adopted distance for the parent galaxies, because 
they depend only on the ratio between the luminosity of 
the echo and that of the SN at maximum. 
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Comparison of the observed spectrum of SN1998bu at phase 670d 
(black line) and with the time integrated echo spectrum computed 
as described in the text (red line). 

The distance D is related to the l inear diamet er of the 
ring through the relation 2R = 2y/ct (2D + ct). In the 
case of SN 1991T the echo has been resolved with a diam- 
eter of 0."4. Based on the Cepheid distance to NGC 4527, 
this results in a linear diameter 2R =120 lyr and hence 
D ~ 150 lyr. Given the observational uncertainties, this 
is in very good agreement with the estimate based on the 
light echo modeling. In the case of SN 1998bu, because of 
the much higher optical depth of the dust, we had to place 
the dust layer at a larger distance than for SN 1991T. In 
turn, this causes the echo ring diameter at any given time 
to be larger. Based on the geometry described earlier we 
obtain that at the time of writing, 2.5 yr after the ex- 
plosion, the diameter of the echo should be ~ 0.3 arcsec. 
Indeed, the echo of SN 1998bu has already been resolved 
using HST (Leibundgut p.c). 

A consistency check of the light echo interpretation for 
the peculiar light curve of SN 1998bu can be obtained 
by comparing the observed and computed spectrum of the 
light echo. Clearly the most important contribution to the 
light echo comes from the light emitted near maximum. 
Spectra of SN 1998bu at eight different epochs from -6 to 
50 days after maximum were published by Hernandez et 
al. (2000). Weighting the spectra according to the inte- 
grated luminosity around each observed epoch and using 
our assumed geometry (cf. Equation 1), we computed the 
contribution of each early-time spectrum to the emerging 
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echo spectrum at phase 670d. The input spectra were cor- 
rected for reddening using a standard extinction law and 
assuming that the scattering function as a similar A -1 de- 
pendence. These corrected spectra were coadded to com- 
pute the expected spectrum of the echo (at 670d the SN 
nebular spectrum is several magnitudes fainter than the 
echo and does not give a measurable contribution). This 
is compared with the observed spectrum in Fig. 3. Given 
the crudeness of some of the assumptions and the incom- 
plete spectral coverage of the near-maximum phase, the 
agreement between the observations and the model is ex- 
cellent, which we consider a strong argument in favor of 
the light echo scenario. 

Moreover, even if we stress that similar results can be 
obtained for different geometries, the fact that with the 
configuration we have chosen it is possible to explain the 
observed light curve, spectrum and reddening is sugges- 
tive. 

3. ARE BRIGHT SN IA LINKED TO DUSTY STAR FORMING 
REGIONS ? 

The observation and interpretation of light echos can 
tell us about dust properties and distribution in galaxies, 
but may have far-reaching applications and consequences. 
One such application is to measure distances of galaxies 
Sparks et al. (1999). 

SN 1998bu is only the second Type la SN for which a 
light echo was observed, the other being SN 1991T in the 
Sbc galaxy NGC 4527. It is remarkable that both SNe are 
very slow declincrs. SN 1991T was also spectroscopically 
peculiar at and before maximum, indicating a higher pho- 
tospheric temperature than in normal SNe la, as shown 
by the presence of strong Felll lines and by the absence of 
Fell lines before maximum (Filippenko et al. 1992). This 
was confirmed by spectroscopic modeling, which also re- 
vealed that SN 199 IT had an abnormally high abundance 
of Fe-group species in the outer, fast-moving part of the 
ejecta (Mazzali, Danziger & Turatto 1995). Such a pe- 
culiar element distribution suggested that SN 199 IT may 
have been the result of an unusual explosion mechanism. 
Further studies revealed that SN 1991T had very broad Fe 
nebular lines Mazzali et al. (1998). This and the bright- 
ness of the SN at maximum led various authors to con- 
clude that SN 1991T produced an unusually large amount 
of 56 Ni, about 1 M (Spyromilio, Meikle, Allen & Graham 
1992; Mazzali, Danziger & Turatto 1995; Cappellaro et al. 
1997; Mazzali et al. 1998) though this may be challenged 
by recent calibration of the SN absolute luminosity (Saha 
et al. 2000; Richtler et al. 2000; Gibson & Stetson 2000). 
Finally, Fisher, Branch, Hatano & Baron (1999) suggested 
that SN 1991T may have been the result of the explosion 
of a super-Chandrasekhar mass progenitor. 

SN 1998bu shares many of the properties of SN 1991T: 
it was bright at maximum, it had a slow decline and broad 



nebular lines, although not quite as extreme in these re- 
spect as SN 1991T. It is therefore very interesting that 
both SNe are significantly reddened and show a dust echo. 
Other two nearby SNe la were heavily reddened: SN 1989B 
was a normal SN la for which indeed there may be some ev- 
idence of a light echo (Milne, The & Leising 1999), while 
SN 1986G was a fast decliner and spectroscopically pe- 
culiar object for which, despite an early claim (Schaefer 
1987), there was no evidence of a light echo (Turatto et al. 
1990). The latter requires that the dust was more distant 
from SN 1986G than in SN 1998bu or SN 199 IT (by plac- 
ing the dust cloud 10 times further the echo magnitude 
becomes 2.5 mag fainter). 

Historically, the fact that SN la were seen as very homo- 
geneous was attributed to their progenitors arising from 
a single stellar population. Adding to this that SN la 
were found in all types of galaxies, even in Ellipticals, con- 
tributed to the standard paradigm that the progenitors of 
all SN la belong to the old stellar population. Now there 
is evidence that slow decliners (or high luminosity SN la) 
may be preferably associated with a younger stellar pop- 
ulation (and therefore conceivably more massive progeni- 
tors), a suggestion first made by van den Bergh & Pazder 
(1992) and confirmed by Hamuy et al. (1996, 2000). Our 
result that both SN 1991T and SN 1998bu are located 
close to dusty regions is coherent with this scenario. 

Obviously one needs to improve the statistics by adding 
more cases. In this respect one interesting object is 
SN 1998es. This is a SN la whose spectrum near max- 
imum was very similar to that of SN 1991T (Jha et al. 
1998) and shows the signature of quite similar interstel- 
lar reddening (Patat et al. in preparation). Therefore 
SN 1998es maybe an interesting candidate to search for 
light echo though, because of the relatively large distance 
(ji = 33.2), the apparent magnitude may result quite faint 
(V~ 24 if the echo has the same intensity as in SN 1991T 
and SN 1998bu). 

4. CONCLUSIONS 

The peculiar late light curve and spectrum of SN 1998bu 
are attributed to the echo from circumstellar dust of the 
light emitted near maximum. Based on a simple modelling 
and on the comparison with SN 199 IT (which also expe- 
rienced a similar phenomcnom) we estimate that the dust 
cloud is located in front of the SN and relatively nearby 
(~ 100pc). Because of the similar echo luminosity but 
much higher optical depth of the dust in SN 1998bu com- 
pared with SN 1991T, we expect that the dust layer is 
more distant and hence the echo ring of SN 1998bu grows 
faster than in SN 1991T. 

The association of dust with SN la, possibly of a spe- 
cific subtype, may have interesting implications for the 
progenitor scenario and prompts for a renewed effort for 
monitoring SN la at late phases. 
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